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Task  1.  Properties  of  Differential  Equation  Models 

In  [18,32,40]  we  studied  stability  properties  of  unstable  fluid  flow.  We  presented  the  effects 
of  the  nonlinearities  and  properties  such  as  bifurcation  in  [18,36,51].  Stability  of  multiphase 
flow  was  studied  in  [32].  This  gives  us  insight  into  the  important  properties  of  unstable, 
bubbly  flow. 

Special  properties  and  treatment  of  singularities  in  the  solutions  to  partial  differential 
equations  were  presented  in  [4.39]  Asymptotic  behavior  of  solutions  was  studied  in  [12,36]. 
Properties  of  hyperbolic  conservation  laws  and  non-hyperbolic  systems  were  analyzed  in 
[17,18,36,37,39,52],  Properties  of  models  of  various  types  were  considered  in  [4,62,63]. 

Finally,  unknown  parameters  must  often  be  determined  in  various  partial  differential 
equation  models.  Special  direct  marching  methods  to  treat  unknown  coefficients  in  time- 
dependent  problems  have  been  developed  in  [31,38,65-68,76]. 

Task  2.  Analysis  of  Numerical  Procedures 

Numerical  properties  of  the  various  discretization  schemes  for  the  partial  differential  equa¬ 
tion  models  must  also  be  well  understood.  Stability  and  accuracy  for  finite  element  methods 
are  studied  in  [2,6,7,30,35,41,56,70,74],  for  finite  strip  methods  in  [19-21],  and  for  colloca¬ 
tion  methods  in  [14,53,54].  Analysis  for  models  for  chemically  reacting  flows  is  presented  in 
[8,15,35,42,72,75]. 

Accuracy  of  fluid  velocities  is  essential  in  the  coupled  systems  in  our  models.  Mixed  finite 
element  methods  have  been  studied  in  [30,34,42,45].  Mixed  methods  have  been  coupled  with 
other  techniques  in  multicomponent  and  multiphase  problems  in  [16,36,40,50,64,81,84,87]. 
The  degeneracies  in  the  multiphase  case  cause  difficulties  in  the  analysis  [16]. 

Hyperbolic,  near  hyperbolic,  and  non-hyperbolic  systems  were  studied  in  [2,5,33],  in  con¬ 
nection  with  convection-diffusion-like  problems  in  [5,54].  Front  tracking  methods  have  been 
analyzed  via  Riemann  problems  in  [36,39,52].  Hyperbolic  equations  are  also  used  in  in¬ 
verse  methods  for  seismic  prospecting,  sonar,  nondestructive  testing,  etc.  Extensive  work  on 
accuracy  of  numerical  methods  for  interface  conditions  was  studied  in  [25,79,89]. 

Operator  splitting  techniques  [2,3,16,72,81,86]  have  been  presented  using  modified  method 
of  characteristics  concepts  [2,8,16,43,73,84].  When  applied  to  multiphase  problems  or  prob¬ 
lems  with  nonlinear  flux  functions,  the  operator  splitting  methods  also  utilize  Petrov-Galerkin 
methods  [2,41,73]  when  the  problems  are  not  fully  symmetrized  by  the  modified  method 
of  characteristics.  Extensions  of  these  methods  to  Eulcrian-Lagrangian  Localized  Adjoint 
Methods  (ELLAM)  conserve  mass  and  treat  boundary  conditions  well  [29, 51,77, 86, SS]. 

Full  asymptotic  error  estimates  are  developed  for  the  difficult  coupled  systems  of  par¬ 
tial  differential  equations  in  [6,13,42,45].  Superconvcrgcnce  results  for  mixed  finite  element 
methods  along  Gauss  lines  were  obtained  in  [34], 

Flow  of  various  fluids  through  heterogeneous  media  can  be  quite  unstable,  with  fingering 
instabilities.  Development  of  effective  parameters  to  model  those  instabilities  for  nonlinear 
scale-up  is  considered  in  [32,78,80,82]. 
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Task  3.  Adaptive  Local  Grid  Refinement 

A  major  effort  has  been  concentrated  on  various  aspects  of  adaptive  local  grid  refine¬ 
ment.  Elementwise  local  refinement  techniques  were  presented  in  [7,53]  for  Burgers’  Equa¬ 
tion.  Adaptive  techniques  that  entail  error  estimators  or  AI  strategies  were  discussed  in 
[1,7,10,16,48,53]. 

Purely  local  refinement  techniques  which  involve  data  structure  trees  have  been  described 
in  [1,10,33].  A  more  efficient  method  that  requires  a  simpler  data  structure  is  a  “patch 
refinement”  technique  [3,40,44,46,48,65,71,82].  Aspects  of  accuracy  of  the  composite  grid 
were  considered  in  [22,23,46,57,61,65]  and  efficiency  of  solution  in  [24,44,46,87].  BEPS  and 
FAC  methods  are  presented  and  compared  in  [9,22,28,46,57,65].  Extensions  of  these  concepts 
to  local  time-stepping  methods  appear  in  [27,60].  Applications  of  these  methods  in  industrial- 
level  simulators  appeared  in  [22-24]. 

The  application  of  patch  method  to  follow  and  resolve  moving  fluid  interfaces  is  presented 
in  [16,41].  Parallel  aspects  appeared  in  [14,55]. 

Several  aspects  of  front  tracking  methods  have  been  studied.  Riemann  problem  analyses 
appeared  in  [36,39,52]  while  the  complications  due  to  non-hyperbolic  systems  were  studied 
in  [18,39,52].  Combinations  with  mixed  finite  element  methods  are  considered  in  [34,45,47], 

Task  4.  Computational  Considerations 

The  advent  of  parallel  and  vector  computer  architectures  and  visualization  capabilities  are 
revolutionizing  computational  considerations.  Basic  efficiency  considerations  were  discussed 
in  [3,10,47,49,59],  Parallel  and  vector  techniques  have  been  presented  in  [1,14,19,21,46,55,83], 
with  scheduling  analysis  in  [1].  Many  of  the  local  grid  refinement  techniques  have  been 
based  upon  domain  decomposition  methods  [3,9,12,16,28,41,44,50,80].  The  potential  for 
parallelization  of  the  modified  method  of  characteristics  was  pointed  out  in  several  papers 
and  is  currently  being  implemented. 

Due  to  the  size  and  poor  conditioning  of  the  systems  involved,  preconditioned  iterative 
methods  are  essential  for  linear  solution.  Preconditioners  are  presented  for  conjugate  gradi¬ 
ents,  conjugate  residuals,  and  domain  decomposition  methods  in  [3,9,44,85]. 

Finally,  we  have  made  substantial  progress  in  visualization  on  our  Ardent  Titan  and  our 
Silicon  Graphics  workstation.  This  greatly  aids  our  computational  capabilities,  allowing  us 
to  use  the  computers  in  an  experimental  mode  to  build  intuition. 

We  emphasize  that  the  combination  of  modeling,  analysis,  numerical  analysis,  and  com¬ 
puting  skills  is  essential  in  addressing  the  complex  and  unstable  coupled  systems  of  partial 
differential  equations  used  to  model  problems  of  interest  to  the  Navy. 
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